Non-Perturbative QED and QCD at Finite Temperature 



Michael Strickland 
Department of Physics, Ohio State University 
OH, 43210, USA 
Email: strickland.4i@0su.edu 

We present results of numerical solutions of Schwinger-Dyson equations for the 
finite-temperature quark and electron propagators. It is shown that both strongly 
coupled QED and QCD undergo a chiral symmetry restoring phase transition as 
the temperature is increased. We go beyond the bare vertex or "rainbow" approx- 
imation by applying the finite-temperature Ward-Takahashi identity to constrain 
the non-perturbative vertex function. 



1 Introduction 

Within both strongly coupled quantum electrodynamics (QED) and quantum 
chromodynamics (QCD) there are analytic and numerical results whiclLindi- 
cate that there is dynamical chiral symmetry breaking in these theoriesEl'El At 
finite temperature it is widely believed that QCD undergoes a phase transition 
from the low-temperature phase with confined quarks and gluons and chiral 
symmetry breaking, to a high-temperature phase in which quarks and gluons 
are deconfined and chiral symmetry is restored. Although the deconfinement 
and chiral phase transitions are two separate transitions, lattice studies find 
that the two transitions occur at approximately the same temperature. Like 
QCD, strongly coupled QED is expected to undergo a finite-temperature phase 
transition which restores chiral symmetry. The main tool for the theoretical in- 
vestigation CLfppon-perturbative finite-temperature QCD has been lattice gauge 
simulations H'H. However, the results from these studies are limited to small 
lattice sizes due to the numerical complexity of lattice simulations. In ad- 
dition, extraction of detailed information (eg, bound-state wavefunctions) is 
complicated by the statistical nature of the data. Another approach is to solve 
truncations of the finite-temperature Schwinger-Dyson (FTSD) hierarchy to 
obtain information about the phase structure of QED, QCD, and other field 
theories. 

In this work, we study solutions of the FTSD equation for the electron and 
quark two-point functions in the imaginary-time formalism. We extend previ- 
ous treatments by using a vertex function which satisfies the finite-temperature 
Ward-Takahashi identity. Our results indicate that both strongly coupled QED 
and QCD undergo a chiral symmetry restoring phase transition. The organiza- 
tion of the paper is as follows: In Section II we present the finite-temperature 
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formalism. In Section III results of the numerical solution of the FTSD equa- 
tion are presented and discussed. Conclusions are drawn in Section IV. 



2 Finite Temperature SD Equations 

In the imaginary time formalism the finite-temperature Schwinger-Dyson equa- 
tion is 

T.{p,u) = S~'^{p,u) - S~'^{p,u) 

-"/-^'->-»^<— 

where = l/^J2mI d'^k/{2-KY, f = g^{N^ - l)/(27Vc), and u is the heat 
bath four-velocity. The sum is over odd Matsubara frequencies.]^ 

2.1 Finite temperature quark and gluon propagators 

When written in covariant form, the most general finite-temperature quark 
propagator is 

S^^{p,u) ^ iA{p^ ,p-u)i) + B{p'^ ,p-u) + iC{p^ ,p-u)i . (2) 

A tensor term proportional to a^u is ruled out by VT invariancei 

The finite-temperature gauge boson propagator in a covariant gauge, as- 
suming q^W-^" = 0, has the general form 

where a is the gauge fixing parameter, P!^^ is the four dimensional longitudinal 
projector, and and Pl^^-wce three dimensionally transverse and longitudi- 
nal projectors, respectivelyt£l 

2.2 Finite Temperature Rainbow Approximation 

If we replace by the bare vertex, in (|^) we get 

Jk 

"The QED FTSD equation is obtained when g — » e. 
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Inserting the general forms of the finite-temperature quark and gluon propaga- 
tors and taking traces of (^) with — and / gives three coupled integral 
equations 

^c^{p^C~{p-u)A)/2\ (5) 
where = A — Z2, Sc ^ B — M^, and = C, and 

A = Tr[ — i'^Y,(jp,p ■ u) ] 
B = Tr[lE(p,p- u)] 

C = Tr[ — I](p,p • u) ] . (6) 
When the traces in are performed we get 

-^=79^ / kdk 



fT ^ p ^ i^B(p,a.„,k,c.„;B) 

^ = 7^ L / kdk 



(2^)2 ^f-^ Jo " -\A^k'+u;l) + 2ACu;,„ +B^ + C^ 
(2^)2 - -\A^k^+ ujI) + 2ACuj,n + B^ + Cy ' 



where Ka, Kb, and Kq are the integration kernels which depend on the form 

1.2 
^3- 



of A3L and A3T and k^ = fcj + A;| + fc| 



^.5 Finite Temperature Quenched Rainbow Approximation 

When the gluon or photon propagator is taken to be_bare (A3L = A3T = 1 = 
A4T) the kernels in equations (0) can be evaluated^ The resulting equation 
can be solved iteratively by evaluating the integrals using Gaussian quadra- 
tures and performing the Matsubara sums numerically. The temperature de- 
pendence of the scalar coefficient function is shown in Figj^. We see that the 
mass is relatively unchanged until T just below the critical temperature for 
chiral symmetry restoration, at which point it undergoes a very sharp transi- 
tion into the symmetric phase. This suggests that the chiral phase transition 
is first order in quenched QED/QCD. 

Note here that the discontinuities in _B as a function of T seen in Fig. || 
are a result of the finiteness of the UV cutoff, A. Since ujn^^^ ~ A, as the tem- 
perature increases fewer modes are included in the Matsubara sum. If the UV 
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Figure 1: B(p = 0, n = 0) as a function of T/A showing the chiral symmetry restoration 
temperature Tc(g^ = 20 , A = 1000). 

cutoff, A, were significantly larger, then these discontinuities would be much 
smaller, and as A oo these discontinuities would be removed completely. We 
are unable to take A very large due to the computer time needed to compute 
the Matsubara sums. 

Using this method, we have calculated the temperature dependence of the 
critical coupling for dynamical mass generation. The result of this calculation 
is given in Fig. |^ showing that the critical coupling increases with temperature. 
Therefore, for fixed coupling, go, above the critical coupling there is a finite- 
temperature chiral phase transition when gc{T) > go- 

2.4 The Finite Temperature Vertex Function 

The requirement that physical quantities be independent of the choice of gauge 
leads to non-perturbative relations between the fermion propagator and the 
vertex fimr-timi \^ Within QCD this relation is called the Slavnov- Taylor 



(ST) identitjtaa 

ip~q),A:ip,q;p^q)[l + biq'')] - 

[gT" - B-{p,q;p- q)]S-\p) - S-\q)[gT- - B-{p,q-p- q)] , (8) 



where a = 1 . . . N"^ — 1 is the color index, b{k'^) is the ghost self-energy, and 
is the ghost-quark scattering kernel. 
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Figure 2: Critical coupling as a function of temperature along with fit of the form gc{T) 
A + BT'^ including the effects of the field renormalizations ( A = 1000) . 



For QED this identity reduces to the abehan version of the ST identities, 
the Ward-Takahashi (WT) identity*" 



i{p- q)vK{p,q;p- q) ^ S ^{p)-S ^{q). 



(9) 



The most general form for the finite-temperature fermion-gauge boson 
vertex function is 

A^{p,k;p-k) = A^{p,k-,p-k)+Al{p,k;p~k) 

= {p + k)t, [a{p, k, u){^ + ^) + b{p, k, u) + c{p, k, 

+d{p, fc, u)-f^ + [e{p, fc, u){-^ + ^) + fip, fc, u) + g{p, fc, u)if\ 
+Al{p,k-p-k), (10) 

Inserting the finite-temperature fermion propagator (||) and vertex function 
( p!o| ) into gives four equations. Unfortunately we have seven unknowns so 
this is not enough to solve for all of the coefficient functions that appear in 
the finite-temperature vertex function. We can, however, find a one-parameter 
solution which satisfies the finite-temperature Ward-Takahashi (FTWT) iden- 
tity: 



a(p, k,u) — a 



A{p'^,p-u) - A{P,k-u) 
2(p2 - k"^) 



''For axial gauges, there are no ghost fields and therefore, within these gauges, the ST and 
WT identities are equivalent. 



5 



b{p, k, u) 


= a 


c{p, k, u) 


= a 


d{p, k, u) 


1 

~ 2 


e{p, k,u) 


= (1 


f{p,k,u) 


= (1 


g{p,k,u) 


= (1 



B{p'^,p-u) - B{k^,k-u) 

i(p'^ — k^) 
C{p ^,p-u)-C{k^,k-u) 

^{A{p',p-u) + A{k',k-u)) 





p-u) - 


A{k^,k- 


u) 




2u-{p 


-k) 




B{p^ 


p-u) - 


B{k'^,k 


u) 




iu- [p 








p-u) - 


C{k^,k 


u) 



-{p~k) 



(11) 



The solutions for e, /, and g are divergent for n = m,p ^ k. Choosing a = 1 so 
that all of the coefficient functions are finite we obtain the following solution 
to the FTWT identity 



a{p, k, u 
b{p, k, u 
c{p, k, u 

d{p, k, u 
e{p, k, u 



A{p'^,p-u) - A{k'^,k-u) 



2(p2- 




B{p^,p-u) - 


B{P,k-u) 


i{p'^ - 


-fc2) 


C{p^,p-u) - 


C{k^,k-u) 






- {A{p'^,p-u 


) + A{k'',k- 



f{p,k,u) = g{p,k,u) = 0. 



(12) 



2.5 Temperature Dependence of the Quark Condensate 

To exactly solve the Schwinger-Dyson equations for the quark propagator we 
must also solve the coupled equations for the gluon propagator and all of the 
higher n-point functions. Short of this we can propose analytic forms for 
the gluon propagator and study the resulting quark propagators. At zero- 
temperature there has been much work along these lines. 

Since perturbative renormalization group studies are not reliable below 
p^ ^ 1 GcV^ the running of the coupling constant at low energy must be deter- 
mined in a different manner. Results from lattice and analytic studies suggest 
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Figure 3: Quark condensate {— <gq^f^ as a function of temperature with the bare vertex 



(o) and the improved vertex (•) Eq. 
10 GeV) 



(y). { (m2 = .69 GeV, g"^ = 17.546, Mo = 1 MeV, A = 



that the zero-temperature gluon propagator is significantly enhanced for small 
spacelike-/c^.^^~^'* Traditional potential analysis suggests that the gluon prop- 
agator behaves like 1 /q'^ for small q and studies have been performed with this 
assumption or regularizations of this form.^^^^^ When modeling confinement 
we will instead use a parameterization which is a regularized version of 
plus a perturbative piece 

— 2 — = 2~, — 2 — ='m't^f^^ [Q)SnO + —, — J, 13 



where d — 12/(33 — 2Nf). This parameterization contains two pieces: a delta 
function which models long-range effects like confinement through an inte- 
grable infrared singularity, and a perturbative piece|-|Which ensures that the 
large-g^ limit is correct up to logarithmic correctionsQ The parameter rrit can 
be interpreted as the energy scale which separates the perturbative and non- 
perturbative regimes. In a zero-temperature calculation using a similar pa- 
rameterization, Roberts and Frank have fixed uit by fitting to a set_pf pion 
observables (eg, iV^ , scattering lengths) and find = 0.69GeVllil 

Once the gluon two-point function is specified, the FTSD equation with 
the vertex function (|o|jl^) gives the temperature dependence of A, B, and C. 
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The quark condensate can be written in terms of these functions 



^- ^_ B{k\LOm) Mo \ 

'^'^ Jk UV + 2ACojm + B^+C^ p2 + ) ■ ^ ' 

In Fig.^ we plot the temperature dependence of the quark condensate 
within the rainbow approximation and with the improved vertex ( p^|p^ ) . The 
magnitude of the condensate and the critical temperature depend sensitively on 
the vertex used, demonstrating that proper treatment of the finite-temperature 
vertex function is necessary. In particular, the rainbow approximation vertex 
function breaks gauge invariance and therefore the critical temperature will 
depend on the choice of gauge. By using a vertex function (|o|jl^) which sat- 
isfies the FTWT identity (|) the gauge dependence of the critical temperature 
is reduced. 



3 Conclusions 

In this work we have shown that is it possible to solve truncated FTSD equation 
beyond the bare vertex approximation. The propagators obtained allowed us 
to study the nature of the chiral phase transition which we find to be first- 
order when the bare gauge boson propagator is used. The coefficient functions 
which appear in the finite-temperature fermion propagator are determined non- 
perturbatively ancLcan be used to calculate other quantities like and 
in the chiral limitu When there is explicit chiral symmetry breaking, these 
propagators can be used in solutions of the finite temperature Bethe-Salpeter 
equation in order to obtain the temperature dependence of the masses and 
widths of the light to intermediate mesons. 
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